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In a recent paper [1] a new initial value formulation of fermionic QFT was 
presented that is applicable to an arbitrary observer in any electromagnetic 
background. This approach suggests a consistent particle interpretation at all 
times, with the concept of 'radar time' used to generalise this interpretation 
to an arbitrarily moving observer. In the present paper we extend this formal- 
ism to allow for gravitational backgrounds. The observer-dependent particle 
interpretation generalises Gibbons' definition [2] to non-stationary spacetimes. 
This allows any observer to be considered, providing a particle interpretation 
that depends only on the observer's motion and the background, not on any 
choice of coordinates or gauge, or on details of their particle detector. Consis- 
tency with known results is demonstrated for the cases of Rindler space and 
deSitter space. Radar time is also considered for an arbitrarily moving observer 
in an arbitrary 1+1 dimensional spacetime, and for a comoving observer in a 
3+1 dimensional FRW universe with arbitrary scale factor a(t). Finite volume 
measurements and their fluctuations are also discussed, allowing one to say 
with definable precision where and when the particles are observed. 



Key Words: particle creation, fermion, observer, Slater determinant, radar time. 



1. INTRODUCTION 

Our recent initial value formulation of fermionic QFT [1] emphasised the 
states of the system, described in terms of Slater determinants of Dirac 
states. The vacuum was defined as the Slater determinant of a basis for 
the span of the negative spectrum of the 'first quantized' Hamiltonian, thus 
providing a concrete manifestation of the Dirac Sea. As well as generating 
simple derivations of the general S-Matrix element and expectation value in 
the theory, the approach suggested a consistent particle interpretation at all 
times. In the present paper we extend this work to encompass gravitational 
backgrounds. 
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In common with the Canonical, Tunnelling, and Wave-functional ap- 
proaches, our analysis depends explicitly on the choice of 'in' and 'out' 
particle/antiparticlc decomposition, and on the corresponding categorisa- 
tion of in and out modes. Various categorisation schemes have been pro- 
posed, based on asymptotic or adiabatic properties of solutions, or on the 
diagonalisation of a suitable Hamiltonian. However, most schemes depend 
on a choice of coordinates or of gauge [3] , and the relation of these choices 
to the motion of an observer or the behaviour of a particle detector is of- 
ten ambiguous [4]. Also, those schemes based on asymptotic or adiabatic 
approximations generate accurate predictions only at asymptotically late 
times or in sufficiently weak backgrounds. Particle detectors provide a more 
operational particle concept, but their predictions are not always propor- 
tional to the number of particles present [4, 5] even when the detector is 
inertial. 

The categorisation scheme proposed in [1] and developed here provides a 
consistent particle interpretation at all times without requiring any asymp- 
totic conditions on the in and out states. It consistently combines the 
conventional 'Bogoliubov coefficient' and 'tunnelling amplitude' methods, 
resolving gauge inconsistencies that trouble each of these [3]. Also, by 
utilising the concept of radar time, the present definition naturally incor- 
porates the motion of the observer (detector) , providing a definition which 
depends only on the observer's motion and on the background, and not on 
the choice of coordinates or gauge. Since it is applicable at all times, we 
can state not only how many particles are created, but also when they were 
created, and how they behaved after their creation (as in the application of 
[1] to spatially uniform electric fields [6]). By considering finite- volume op- 
erators with controllable fluctuations, we can further specify, with definable 
precision, where the particles are created. 

In Section 2 we specify the representation of states used, and their evo- 
lution, and present formulae for the general S-Matrix element and the ex- 
pectation value of the theory. The particle interpretation is described in 
Section 3, specifying the states of Section 2 in terms of their particle con- 
tent. Finite volume measurements are also discussed in this Section, along 
with their fluctuations. For concreteness we suppose that the spacetime is 
globally hyperbolic, and that the observer-dependent foliation is Cauchy. 
Strictly, this excludes situations with particle horizons, but we present a 
simple example in Section 4 which shows that the formalism is still well 
suited to the treatment of horizons. We consider massive and massless 
Dirac fermions in 1+1 dimensional flat space as seen by a uniformly ac- 
celarating observer, and demonstrate consistency with known techniques, 
by rederiving the well-known thermal distribution of Rindlcr particles. We 
consider the spatial distribution of these Rindler particles, and comment 
briefly on fluctuations. 
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In Section 5 we present examples of radar time, treating an arbitrary 
observer in an arbitrary 1+1 dimensional spacetime, and also a comoving 
observer in an FRW universe with arbitrary scale factor a(t). DeSitter 
space and the Milne universe are described in more detail. We conclude 
with a brief discussion in Section 6. An appendix describes the connection 
between projection operators and 2-point functions, and emphasises the 
role of the negative energy Wightman function as the "Dirac density matrix 
of the Dirac Sea" . 

2. THE STATE SPACE 

2.1. Preliminaries 

The Lagrangian density for the Dirac equation in gravitational and elec- 
tromagnetic backgrounds is [7, 8] 

C = $[e{x)i;{x)(i^(x)V^ - m)V(ar)] (1) 

where 7 M (x) = e£(x)j a satisfies {"/^(x),^ (x)} — 2g tJ - u (x), j a are a repre- 
sentation of the normal flat space Dirac matrices, i/j = ^7°, e%(x) is the 
Vierbein, and e(x) = det(e^(a;)) = ^J— det(g^ v (x)). The covariant deriva- 
tive Vft acts on spinors according to: 

\7 tl ,p(x)^(d fl + T ( , + ieA ( ,)iP{x) 

where T M = \"fv^ ^f u , is the electromagnetic potential, and e is the 
charge of the fcrmion (e<0 for electrons). is often written [9] in terms 
of the 'connection field' uj^ b (x) as = — jui" <J a b where a a b = \ [la , lb] • 
This Lagrangian gives rise to the governing equation: 

(i^{x)\7 f ,-m)^(x) = (2) 

and the energy-momentum tensor: 

T^)=^ 1{l y v) ^]-^f\c (3) 

eyx) 

The inner product (V , | < ^)e on the spacelike Cauchy surface S is given by: 
(V#)e = J e(x)TP(x)^(x)<f>(x)d^ (4) 
and is independent of £ by virtue of (2). 
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The 'first quantized' state space on some spacelike hypersurfacc 

£ can be defined as the space of all (finite norm) spinor valued functions 
of x|s (the projection of x* 1 onto the hypersurface £). We restrict our 
attention to hypersurfaces £ that are Cauchy, so that the various 7i(S) 
are all unitarily equivalent, and can simply be denoted H. We will denote 
a first quantized state on £ by tp(x\^) or \ipT.) or, where no ambiguity is 
possible, simply ip. (There will be little need to distinguish between a state 
and its coordinate representation.) 

2.2. The Full Fock Space Over H 

The antisymmetric Fock Hilbert space over the complex Hilbcrt space 
H is denoted T^{TL) and is defined [10] in terms of the antisymmetric 
Tensor Algebra over H. It is a natural and familiar construction by which 
a quantum theory of fermions can be formulated. This construction is 
described in [1], and as a preliminary we outline it here. Let H be the 
Hilbert space in the previous Section, with inner product denoted by ( | ). 
Let ® n TL denote the direct product of n copies of H, and let A n H denote 
the restriction of ® n 'H to those states which are completely antisymmetric 
under changes in the order of the elements <G H from which the state 
is constructed. Given \ipi) \tp 2 ) ■ ■ ■ G ® n W, we can define A V>2 A • • • A 
%l> n ) € A n H by: 

I V>i A V>2 A • • • A %j) n ) = -p= V] sign(er) | ) |VV(2) ) ■ ■ ■ I VV(™) ) (5) 

where {a(i),i = l,...,n} is a permutation of {l...n}. This is simply the 
Slater determinant of the states \ipi) . . . \ipn)- The antisymmetric Fock 
Hilbert space is now given by: 

where A°H = C and A 1 H = H. States which lie entirely within h r 7i for 
some r are said to be of grade r. 

A useful operation on !F/\(H) is the 'inner derivative' i^ : A n H^> /\ n ~ 1 H 
(named by analogy with differential geometry). This is defined by: 

i^\(j) l A---A^ n ) = ^{-) l+1 {i>\^)\(j} l A---A^A---A^ n ) (6) 

i 

where the check over 4>i signifies that this state is omitted from the prod- 
uct. The relation i^ : JF A (7i) — > J 7 A (Tl) is obtained from (6) by imposing 
linearity, together with the convention i^\ = for Ag A°H. It is clear that 
i^(i^\F)) = for all \F) e T A (H), and that: 



iip(cj>A\F)) = ^\<P)\F)-<pA{H\F)) 



(7) 
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The operation i$ is denoted a(tp) by Ottesen [10], and plays the role of an 
annihilation operator. Here i^ will play a similar, although not identical 
role. 

Finally, the inner product on T A (7i) is given by: 

(^1 A • • • Aip n \<t>i A • • ■ A(j) m ) = 5„ m det[(^i|^>] (8) 

where (ipi\(f>j) refers to the inner product on 7i. (For states A,/i <G A°7i define 
(A|/u) = Xfi and (X\F n ) = for any state \F n ) of grade n > 0.) This definition 
agrees with the inner product defined in terms of Slater determinants. 
Although we use the notation ( | ) to refer to both the inner product 
on TL and the inner product on J r A (H), it will always be clear from the 
context which is involved. 

2.3. Operators on Fock Space 

Let A\ :H^H be an operator on the space of Dirac states. We wish to 
construct from it an operator which can act on all of state space. There are 
two useful ways of doing this: Hermitian extension Ah ■ ^a(W) — > -F A (W), 
and Unitary extension Ajj : T A (H) — » -F A (W) (outlined also in Ottesen [10] ) . 
These are respectively defined by: 



A H I ipi A tp 2 A • • • A ip n ) ee \i>i A • • • (ii V><) A ip i+ i ■ ■ ■ A ip n ) (9) 
ic/IV-i A ifa A- • • A Vn) = l(iiV'i) A (A1V2) A • • • A (iiVn)} (10) 

If Ai is (anti)hermitian with respect to the inner product (4) on H, then 
Ah is (anti)hermitian with respect to the inner product (8) on J : A (H). If 
A\ is unitary, then so is Ajj- Also (e Al )u ~e A " , so that if U\ — e Al on 7i 
then U v = e AH on .F A (W). 

Some Simple Properties 

1. (1 + 5)^ = 1^ + ^, [i ff ,B H ] = [i,B] ff and (AB^ = i c/J B [/ - 

2. [i H ,'i/ , A] = (iiV') and [A H ,i^] = -i A ^ 

3. If , tp2 , ■ • • V'n are all eigenstates of Ai with eigenvalues Ai , . . . A„ , then 
|V>i A-02 • • • A ipn) is an eigenstate of Ah with eigenvalue 
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4. If ipi,ip2, ■■■fpn are orthonormal and \F) = A ^2 • • • AV>n) then 

n 

(F\A H \F)=J2mAi\i>i) (11) 
i=l 

n 

^|(i ff ) 2 |^)=^(^|i?|^)+2^((Vi|ii|V'i)(V'i|ii|^) (12) 

-(^|ii|^)(^-|ii|^)) 
(F|(i ff ) 2 |F)-«F|A ff |F)) 2 = X;^il^l^>-El^l^l^>l 2 ( 13 ) 



2.4. Evolution of States 

Given two Cauchy surfaces Si and S , we can define the evolution op- 
erator J7i(Si,S ) on H by 

^i(E 1 ,E )|^ > = |^ Eo (Si)> (14) 

where \ipx ) represents some chosen initial conditions i[)(x\e ) on S , and 
IV'So (^1)) represents the corresponding solution ip(x) of the Dirac equation, 
expressed on the hypersurface Si. 

We consider only QFT in an (external) gravitational or electromagnetic 
background, so that we ignore direct particle-particle interactions and work 
within the 'zeroth order Hartree-Fock' approximation. This assumes that 
the evolution operator on T^iji.) is just the unitary extension of the evo- 
lution operator on H. The action of J7(Si,S ) is now given by: 



f/(Si,So)|Vi,E A • • ■ A Vn,E ) = |V>i,£ (Si) A • • • A V„,e (Si)) (15) 

The multiparticle solution is simply the Slater determinant of the appro- 
priate 'first quantized' solutions. This construction preserves grade, and 
implies that the unitarity of J7(Si,S ) follows immediately from the uni- 
tarity of the first quantized Dirac equation. 

We have now set up a state space, an evolution equation (15) and a 
conserved inner product (8). These are all we need to calculate arbitrary 
S-Matrix elements (from (15) and (8)), arbitrary expectation values (from 
(11)), and even fluctuations in these expectation values (from (13)). How- 
ever, the theory is not invested with physical meaning until the states of the 
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system can be specified in terms of their physical properties. For this pur- 
pose a particle interpretation is invaluable. We now introduce an observer- 
dependent particle interpretation. This will rely on the observer-dependent 
foliation of spacetime provided by Bondi's Radar Time [11, 12, 6]. 

3. OBSERVER DEPENDENT PARTICLE 
INTERPRETATION 

3.1. Bondi's Radar Time 

Consider an observer travelling on path j-.x^ = x ij '(t) with proper time 
r, and define: 

t + (x) = (earliest possible) proper time at which a null geodesic leaving 

point x could intercept 7. 
t~ (x) = (latest possible) proper time at which a null geodesic could 
leave 7, and still reach point x. 
t(x) = ^(t + (x) + t~ (x)) = 'radar time'. 

p(x) = h(r (x) — t~ (x)) = 'radar distance'. 

E To = {x : t(x) = tq} = observer's 'hypersurface of simultaneity at time To'. 




FIG. 1. Schematic of the definition of 'radar time' t(x). 

This is a simple generalisation of the definition made popular by Bondi 
in his work on special relativity and A;-calculus [11, 12, 13]. By construction 
radar time is independent of the choice of coordinates (since no coordinates 
need be introduced in defining it) and it depends only on the motion of 
the observer. It agrees with proper time on the observer's path, and is 
invariant under 'time-reversal' - that is, under reversal of the sign of the 
observer's proper time. 
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It is clear from the definition of radar time that S Tl lies to the future of 
E TQ (for n > ro) except at the observer's particle horizon (if one exists), on 
which the various S T converge. When horizons are present, the domain of 
t(x) is no longer all of spacetime - only that part with which the observer 
can send and receive signals. 

Define now the 'time-translation' vector field: 

w ^jr^ as) 

dx" dx" 

This represents the perpendicular distance between neighbouring hyper- 
surfaces of simultaneity, since it is normal to these hypersurfaces and it 
satisfies k^{x)-§^ = 1. Now use the identity iff^V IJ/ = ik^V IJ/ + a^ v k IJ ,V l/ 
(where k = fc M 7 M and a^ v = |[7 P ',7"] = a ab e%el) to write the Dirac equation 
as: 

^fc^V A^ ^ = H nh (T)^ = -a^k^^ + mk^ (17) 

Here H n h(T~) is not in general Hcrmitian! (hence the subscript nh). At 
first sight this seems to disagree with unitarity on TL. However there is 
no inconsistency, because the inner product now depends explicitly on r 
(via the volume element on E r ), and since (17) is no longer of the form 
i4:\ip(i)) = Hi(t)\tp(t)) , invalidating the standard equivalence proof of uni- 
tary evolution and a Hermitian Hamiltonian. 

To investigate the relation between H n h(To) and the energy momentum 
tensor, define: 

H T M)= f T^(x))k»d^ (18) 
Substitution of (3) into this gives: 

F ro (V) = 3i[(V|^„h(ro)|V)E TO ] = (V|ffi(To)|V)E TO (19) 
where H^tq) = \{H^) +h\{t q )} 

We can define the projection operators : TL — > (to) , and hence the 
spaces TL ± (to), by requiring that are orthogonal projections satisfying: 

H T0 {P+J)> H T0 (V) > H T0 (P-J) (20) 
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for all ip£H, as in [14]. This definition depends only on the background 
and the motion of the observer, and not on the choice of coordinates or 
gauge. It is equivalent to defining: 

H + (to) as the span of the positive spectrum of Hi{t ) 
H~(to) as the span of the negative spectrum of #i(to) 

Here 7i + (ro) is the set of all positive energy states, and H~(to) is the 
set of all negative energy states as defined on S T0 . This definition gen- 
eralises Gibbons' approach [2] to arbitrary observers and non-stationary 
spacetimes. 

In the Canonical approach to background QFT, this split of Ti into posi- 
tive/negative energy modes can be achieved by Hamiltonian diagonalisation 
of the second quantized Hamiltonian that arises by substituting the field 
operator ip(x) into expression (18) for H To {ip). Hamiltonian diagonalisation 
has been criticised [15] for its reliance on the choice of a second quantized 
Hamiltonian which depends on an apparently arbitrary choice of hypersur- 
face and time-translation vector field k^ . Here this arbitrariness has been 
resolved (and in [1, 16]) by specifying the hypersurface S ro and the vector 
field in terms of the worldline of the observer (or particle detector). 

We have assumed here that H\ (r) has no zero energy eigenstates (for any 
t). Even for inertial observers, there exist topologically non-trivial back- 
grounds for which zero energy eigenstates exist, leading to the existence 
of fractional charge; see e.g. [17]. Although such situations are straight- 
forward to describe within the present approach, we will not discuss them 
further here. When there is a particle horizon, so that the observer's hy- 
persurfaces are not Cauchy, zero-energy eigenstates are plentiful and cor- 
respond to states that are unobservable by that observer. For concreteness 
we will suppose in this Section that all S T are Cauchy. In Section 4 wc 
demonstrate with a simple example that this approach remains useful in 
the presence of horizons. 

3.2. Particle States and S-Matrix Elements 

Having defined H ± (t ) for any given observer, we can now define their 
vacuum |vac ro ) on S To , to be the Slater determinant of any basis of H~ (to) 
(normalised so that (vac To |vac To ) = 1). This specifies |vac ro ) up to an arbi- 
trary phase factor. It is the state in which all negative energy degrees of 
freedom are full, and hence is a concrete manifestation of the Dirac Sea. 

To illustrate this, let {ui yTQ ;i G /}, {vi tTo ;i G /} be orthonormal bases for 
7Y + (t ) and H~(to) respectively, where / is some countable index set (the 
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uncountable case introduces no complications). The vacuum on S TQ can 
be written as: 

|vac To ) = |wi iTo A-y 2 ,ro A...) (21) 

and is independent of the choice of basis for H~(t ) (up to a phase factor) 
because of the complete antisymmetry of the Slater determinant. The 
vacuum on E Tl at some 'time' t\ can similarly be written as: 

|vac tl ) = |ui lTl Av 2 , Tl A...) 

The observer who prepares the state 'at time To', on his hypersurface of 
simultaneity S To , needn't be the same observer who measures that state 
'at time n' (on her hypersurface S Tl ). A typical case in which these two 
observers differ is the Unruh effect, where |vac To ) is the Minkowski vacuum, 
prepared by an inertial observer, and |vac Tl ) is the Rindler vacuum, defined 
in the frame of a uniformly accelerating observer. 

The evolved state |vac ro (ri)) obtained by evolving |vac To ) from E To to 
S Tl is (since U = U\,u)'- 

|vac To (ri)) = |ui )To (Ti)Au 2 ,To(Ti)A...) (22) 

where Vi,T ( T i) denotes the state (called |^s T0 (S Tl )) in (14)) obtained from 
v i,T by evolution to £ Tl . It will not in general be contained in W _ (ti). We 
will often refer to |vac ro (ri)) as the 'evolved vacuum', although it is not in 
general a vacuum state. 

From (8), the vacuum- vacuum S-matrix element is simply: 

(vac ri |vac To (n)) = det[(i>j ;Tl \vj. To 

(n)>] (23) 

The probability that |vac To (ti)) will be vacuum at time t\ is then 'Pyac-^ac = 
|(vac Tl |vac To (ri))| 2 . Although this result can be derived by a number of 
methods once a particle interpretation is specified [18, 19, 20], we believe 
that the present derivation (and in [1, 16]) is clearer and more economical. 

States containing particles can be treated just as easily. For instance, 
a one-electron state (at time To) is of the form u TQ A |vac TQ ), and a one- 
positron state (at time To) is of the form i Vt |vac To ). As expected, electrons 
arc represented by the presence of positive energy degrees of freedom, and 



FERMIONIC QFT IN GRAVITATIONAL BACKGROUNDS 11 



positrons by the absence of negative energy degrees of freedom (note that 
v To A |vac To ) = = i Uro |vac ro ) for all u T() e H + (t ) and v TQ E H~(t )). We 
introduce the symbol K}^"'^) ) to denote an 'in' state of m particles 
(in states Ui t . . - Ui m with i\ < 12 < ■ ■ - i m by convention), and n antiparticlcs 
(corresponding to the absence of states Vj 1 ...Vj n ), prepared at time r . 
This state is given by: 



^hh'-fJr) = H^K.rc A • ..«i m ,7t, AUl.rc A • ..V h , T0 ■ ■ -Av jn , T0 ■ ■ ■) (24) 



where the check over Vj jT0 signifies that this degree of freedom is missing 
from the state, and J — ^ {n+l) + Y^!k—i3k appear as an unimportant sign 
convention. 

The general S-matrix element can immediately be written as: 



X ^3 1 3 2 ---] nl >^ ^3l32---3n> T „ y LJ/ 



= H 



J- J' 



det 



7i* 



m' 



\missingJ \missinp/ 



missing/ 



An - 
%,i •• 

en 

V-m isstnqJ 



JX-3n \ 
\missingJ 



( 3i-3n 1 
\missing, 
. / h—3n \~ 
\missing) 



(25) 



(if m — n = m,' — n', and zero otherwise), where 



aij(Ti,T ) = (u i>Tl \u j>To (T 1 )) 

Pij(n,T ) = {Ui, Tl \v jtTo (T 1 )) 



Hj (n , t ) = (wi, ri |uj )T0 (n)) (26) 

£y (7-l,To) = K, Ti Kto(ti)) (27) 



are the time- dependent Bogoliubov coefficients. The Bogoliubov conditions 
follow from unitarity of the 'first quantized' evolution matrix 



£1(7-1,7-0): 



a(ri,r ) /3(ti,t ) 
7(ti, t ) e(ri,r ) 



The ease with which the general S-Matrix element (25) has been derived 
contrasts with many conventional formulations [21, 22, 20, 23], as discussed 
in detail in [1] . The reason is in the concrete representation of states given 
by equations such as (21) and (24), and the simple evolution equation which 
allows us to deduce equations such as (22). In Canonical approaches to 
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QFT in a classical background [7, 33, 23] the states are denned implicitly, 
by the requirement ai|vac) — — b;|vac), where the creation/annihilation 
operators aj,i>j are defined implicitly by the CAR's. The derivation of 
S-Matrix elements then involves more round about methods. One such 
method, analogous to that used in [23], is described in [1] and contrasted 
with the derivation above. 

3.3. Expectation Values 

Given an operator A\(t) :H^>H, we can define its physical extension 
A vhys (T):T A {H)^T A {H) by: 

iphys (t) = i ff (t) - (vac T | A H (t) | vac r ) 1 (28) 

This is the relativistic equivalent of the 'one-particle operator' of mul- 
tiparticle quantum mechanics [30], and is expressible as a normal-ordered 
bilinear of the field operator ip(x). This vacuum subtraction is equivalent 
to normal ordering with respect to the particle interpretation at the time of 
measurement. This choice is also made in previous 'Hamiltonian diagonali- 
sation' procedures [24, 25, 26, 27, 28], and uniquely guarantees the positive 
definitencss of H p ^ ys (T) while maintaining (vac T |i/ p h ys (T)|vac T ) =0. This 
is discussed in more detail in [1], where the relation of our approach to 
Hamiltonian Diagonalisation is also described. 

The expectation value of A p h ys (T) in the physical vacuum |vac T ) at time 
r is zero by construction. Its expectation value in the 'evolved vacuum' 
|vac To (Y)) is in general non-zero, and takes the form 

(vac TQ (n ) | iphys (n ) I vac To (n )) (29) 

N N 

= £K7 tt (n)|i4i(Ti)|u i , 7t) (Ti))-5^<w i , ri |i4i(Ti)|i; < , T1 ) (30) 

i=l i=l 

= Traced A++ - 77 f A +e$A+- + /3e 1 ' A~+) (31) 
where we have defined: 

A tk= («j,ril^i(Ti)K,n) ^7fe~ = Kn|ii(ri)|?;fe, Tl ) 

^" = ( w j>i|ii(n)|wfc,n) and AT+^^ji^n)!^) (32) 



= Aj^j if A\ is Hcrmitian 

The relation (ee^)kj —Skj = — {l^)kj has been used in deriving (31) 
from (30). This step relies on the fact that we are vacuum subtracting 
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with respect to the vacuum at the time of measurement. Notice that if A\ 
is conserved at the level of the Dirac equation, then: 

(vac To (n) | A phys |vac To (n)) = Tracc(A (r ) A (n)) 

which can be non-zero even when Ai is independent of time, because of the 
varying particle interpretation. Herein lies an elegant physical description 
of quantum anomalies, as outlined in [1]. A treatment of the axial anomaly, 
described in terms of this physical mechanism, has been given in [29] and 
in [17]. 

The derivation of (F To (ri)|y4phys(Ti)|-PVo( T i)) f° r an arbitrary state |-FV ( T i)) 
is identical to the derivation of (31), with result: 

( (i^- "±) ^ (-1 ) I ^Phys (-! ) | (}^V±) ^ )> 

m n 

(n ) I M (n ) | u ik >ro (ti ) ) - {v jk iT0 (ti ) | A-t (ri ) I v jk >ro (ri ) ) 

+ (vac T0 (n )|i p hy S (Ti ) | vac T0 (n)> (33) 
We now consider some simple examples of physical extension. 

The charge operator is easily recognised as the physical extension of 
the unit operator on Ti, Q = el p hy S - This is appropriate, since the norm 
(ip\li\il)) of a state tp e H is the Noether charge of the Dirac Lagrangian 
that is conjugate to changes in phase. Charge conservation now follows 
directly from the fact that evolution is grade-preserving. 

The operator -/V p h ys (T) that represents the number of particles (includ- 
ing antiparticles) present at time r is the physical extension of Ni(t) = 
-P+ — P~ , where : H — > (r) are the projection operators onto H (r) . 
Clearly Ni(t ) commutes with Hi(to), but does not in general commute 
with time evolution (since it does not commute with H\{t) for t^t ). 
-^phys(To) inherits both of these properties. Therefore the number operator 
Aphys^o) represents a well-defined physical observable, but one which is 
not conserved. When acting on states in standard form, it gives: 

^-^)!Citt),H("'+'0!C:tl),,) ( 34 ) 

so that it is positive definite and has integer eigenvalues. From (31), the 
expectation value of iVph ys (Ti) in the 'evolved vacuum' |vac ro (ri)) is given 
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by: 

(vac To (n ) | iV phys (n ) | vac To (n )) 
= Trace(/3/3 t +77 t ) (35) 

= E^o(n)+^ (n)} 

i 

u i,T% \ v j .Taiji))^ is the expectation value of 
the physical extension of \iii iTl ) (iii iTl |, and represents the probability that 
the degree of freedom u^ Tl is occupied in |vac To (ri)), i.e. that particle i is 
present. A^ To (ti) = (77 f k = J2j \( v i,ri K> (n))| 2 is the expectation value 
of the physical extension of — \v itTl ){v itTl |; it represents the probability 
that the degree of freedom Vi iTl is unoccupied in |vac To (ri)), i.e. that 
antiparticle i is present. The Bogoliubov conditions imply Tracc(/3/3^) = 
Trace(77^), which again expresses charge conservation. By defining the 
projection operators P^(t) by: 

(r) = J2Kro(r)) (Ui,r (r) | - ft (r, ro)P+ u\ (r, r ) (36) 

i 

P-(t) = J2 KM) (Vi,r (r)| = ^(r.TBjP-^^Ttt) (37) 
we can write Nf - (n) as: 

= K.rJ^Cn)!^) ^ T0 (n) - (^, ri |P+(r 1 )|^, ri ) (38) 

We can rewrite (31) in terms of P^(t) as: 

( vac T0 (ti ) | A phys (n ) | vac ro (n )) = Trace( Ax (n ) (P- (n ) - P" ) ) (39) 

= -Trace(A 1 (r 1 )(P+(r 1 )-P+)) 

where we have used P+ (ti ) + P~ (ri ) = 1 = P+ + P~ . This result can be gen- 
eralised to any state \F) of the form |P(r)) = |V>i(t)AV>2(t) • • -A^ n (r)) where 
^i( T ), 1 /'2(T),...V'n(''") are orthonormal (which includes all | ( nl2 '" 4 " 1 ) (r))). 
In this case, define: 

5>i(r)Wi(T)| 



and P+ ) (t) = 
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where ^2 ieI runs over {ipi, . ..ipn}, and J2i<£i runs over the orthogonal com- 
plement of this. Then (33) can be written as: 



(F(r)| J 4 phys (r)|^(r))=Tracc(A 1 (r)(Pp F) (r)-P;)) 

= -Trace(A 1 (r)(P+ ) (r)-P+)) (40) 

|F)( T ^ t |F)' 



= Trace(A 1 (r) (P™ (r)P+ - P+ (t)P" )) 



The relation of these projection operators to the various Greens functions 
of the theory, and to the first order 'Dirac density matrix' of multiparticle 
quantum mechanics (see for instance [30], pp 9-10) is described in the 
Appendix. 

3.4. Fluctuations 

Fluctuations in expectation values can be calculated using (13). In the 
notation of the previous subsection, we can write (13 as: 



(^(r)|(i phys ) 2 |F(r))-((^(r)|i phys |^(r))) 2 ^ ]T \{Ur)\M^{r))t 
= Traced " , (r)ir (t)P+ } (r)i x (r) ) = - lTracc( [P" ) (r) , A, (r)] 2 ) 
In the evolved vacuum |vac ro (r)) this becomes: 



(A 2 ) = (vac To (r)|(ip hys ) 2 |vac To (r)) - (vac To (T)|i phys |vac T (r)> 2 
= ^|K T0 (r)|i 1 (r)|^ T0 (r))| 2 

i,j 

= Trace(P-(r)i 1 (r)P+(r)i 1 (r)) 

= Tracc(A F A F ) (41) 
where A F = ^ A ++ a + e t A~ 7 + /3 f A + ~j + J A + a (42) 

Consider for example fluctuations in N + (t) , the total number of particles 
(not including antiparticles) in the evolved vacuum. Now (r) = P+ so 
that Ap = f^OL, and: 



{{N+f) = Traced (1-/3/3*)) 



(43) 
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Consider a spatially uniform case, with /%>A;qo-(T,To) ~ /3pS\ a (2Tr) 3 S(p— q). 
Then: 



N ^c,ro( T ) _ f d 3 p , 



V 

2 



f d p \B I 



(44) 



while <(%^) > = |/(^W(1-I^l 2 ) ( 45 ) 

This result suggests that for large volumes, and for |/3p| 2 << 1, the fluc- 
tuations in the average particle density in a volume V vary inversely with 
V, and that in a given volume the rms error in the total number of particles 
is proportional to the square root of the number of particles, as would be 
expected for a weakly interacting system. The same result holds for the 
total number of antiparticles, or for the total number of pairs. (The fluc- 
tuations in the total charge are zero, due to the Bogoliubov conditions.) 
However, 'V in (44) and (45) is formally infinite; V = (2n) 3 5 3 (0). Before 
we can confirm our interpretation of (44) and (45) we must consider finite 
volume measurements. 

3.5. Finite Volume Measurements 

To measure the quantity A in a volume V (on E T ) it seems reasonable to 
consider the operator Ay iP hy S obtained by Hermitian extension and vacuum 
subtraction from the operator: 



lv,i = !( 



rA\ + Ai^y)where 9vip(x T ) : - 



ip(x T ) x T eV 

X T iV 



(46) 



where x T is shorthand for x\^ T (the restriction of x onto S T ). The matrix 
elements of Ay,i are related to those of A\ by restricting the volume inte- 
gral to V, taking the appropriate combination of surface terms to ensure 
Hermiticity. In the Canonical approach, Ay.H is similarly obtained by re- 
stricting the bilinear J E e(x)^[Ai^(x)j fl (x)^(x) + i^(x)j ,1 (x)Ai^p(x)]dT lfl 
to the volume V (here ip(x) is the field operator). 

As defined, however, Ay^hys gives rise to some serious problems: 

1. The fluctuations in N+ phys (i-), N~p hys (r), Ov,phys and Hy phys (T) are 
infinite for any finite V, even in the physical vacuum |vac T ), and even for 
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an inertial observer in flat empty space! To see this, note that: 

(vac T \(N± phy f\vac T )=J2\(u j!T \v i , T ) v \ 2 (47) 

ij 

(vac T |(i? y , phys ) 2 |vac r ) = Y,{Ef - Erf\{uj tT \v itT ) v \ 2 (48) 

where (uj. T \vi^ T )v signifies that the integral in the inner product has been 
restricted to V (it would be zero otherwise). For convenience \uj tT ) is 
chosen to be an eigenstate with eigenvalue Ef and \v ijT ) is an eigenstate 
with eigenvalue — . For an inertial observer in flat 1+1 dimensional 
space, and V being the region \x\ < ^, this becomes: 



(vac T |(i?y iPh ys) 2 |vac T ) = J 7^7^ F (P><l) 



wWP( „„ HF 2 (E p +p)(E q + q) 2 sm 2 [(p + g)L] 

^eF(p, q ) = (E p -E q ) _ (I- {Ep+p)[Eq + q) ) {p + q)2 

which is easily seen to diverge (the divergence is just as bad in 3+1 dimen- 
sions). 

2. The expectation values of Ny h (t) in the evolved vacuum |vac T0 (T)) 
are given by: 

(vac T0 (T)|N+ phys (T)|va^ 

(49) 

(vac T jT)|Nv !phys (T)|vac T0 (T))=Tracc(77 t (v iiT |v j!r )v-3fJ(e^ t (u i!r |v j ^)v)) 

(50) 

The first terms represent the obvious contribution from the sum over cre- 
ated particles of the probability distribution of each created particle. The 
final term has no obvious interpretation. It is comparable in magnitude to 
the first term, but is not positive definite, and it affects (vac ro (r) |Ny phys (r) |vac T 
with opposite sign to (vac T0 (r)|N~ phys (<r)|vac T0 (T)). 

All these problems stem from the fact that 0y does not commute with 
. This is the well-known fact that particle states cannot be confined to 
a finite volume without introducing a negative-energy component. Conse- 
quently, even if [Ai,P^] =0 the same will not be true of Ay,i- Hence, even 
if |vac T ) is an eigenstate of A p h ys , it will not in general be an eigenstate of 
j4y iP h ys for finite V. This can be seen in the Canonical approach by consid- 
ering the transition from operators expressed in terms of ip{x) to operators 
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expressed in terms of creation/anihilation operators. The absence of terms 
proportional to aj&j relies on the relation (ui\Ai\vj) =0. However, when 
the integral is restricted to a finite volume this matrix element can be non- 
zero even if \vj) is an eigenstate of A\. It is these finite volume overlaps 
between positive and negative energy states that contribute to (47) and 
(48), and lead to the unphysical second term in (49) and (50). 
This problem is easily overcome. Define: 

% = P+9 V P+ + P-9 V P- (51) 
and = 5(^1 +^i^v) (52) 

The expressions considered in (46) - (50) will henceforth be denoted Ay^} vc , 
^v?phys> etc - A v ,i now has the property that [A Vtl ,P±] =0 (for all V) if 
and only if = 0. Hence, if |vac r ) is an eigenstate of A p h ys , then it 

is an eigenstate of Ay^hys for all V, so that the fluctuations of A in any 
finite volume in the physical vacuum |vac T ) will be zero. This applies for 
instance to N ± (t), Q and H, and resolves problem 1 above. 
If [ii,P T ± ]=0, then: 



(vac To (Ti)|iy,p hys (Ti)|vac To (ri)) = Trace(/3/3 f — A v ) 



where A++ k = (u^A^^Uk^) and Aj~ = (v^AA^in)^). 
The terms describing finite- volume overlaps between positive and negative 
energy states have been removed, resolving problem 2. For instance, we 
can now write: 



^ ;T0 W = ^(<^ > r|^P-P+WI^ > r)=^5R[Tr a ce((77 t )ik(v k , r |v i , T )v)] 

k 

^. W-^(("^|^P+^(T)| U ,, T )-^3?[Tracc((/3/3 t ) ik K T |u i , T )v)] 

k 
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from which it follows that: 

Ny (t) -Tracc( 7 7 t (v k , T |v i:T )v) 



where 
and 



= / e{x)J-»{x)d^ = f n- (x)d 3 x 
Jv Jv 



(53) 



d t 

n T a ( x ) = e ( x r)J Ta l "( x T)-^:\x=x T where x£ = (t T (x),x) 
N Vt ( t )= Traced (u k , T \Ui,r)v) 



n\ ' (x)d 3 x 



(54) 



where 
and 



dnT, 



n+(x) = e(x T )J+"(x T )-^\ x= 



Then Ny T0 (T) = J v e(x)J+'»(x)+J->»(x)d^ and Q VtV) (T) = ef v e(x)J+"(x)- 
J~ ' fl (x)d'E l _ l . The fluctuations in Ny To (r) are given by: 

((iV+) 2 )=Tracc((/3/3 t ) ij (u jiT |u kiT )v(l-/3/3 t )ki(u 1 , T |u miT )v) (55) 

which clearly reduces to (43) when V covers all of S T . In the spatially 
uniform case, with (3px-q,j(T,To) = (3pSx cr (2Tr) 3 S(p—q), it is straightforward 
to show that (45) also holds approximately for large finite V , confirming 
the interpretation of (45) given at the end of Section 3.4. Similarly the 
fluctuations in Ny TQ {r) are given by: 



((N v f) =Tracc((77 t )i j (v j!r |v k ^ T )v(l-77 t )ki(vi, T |v m;T )v) 



(56) 



In the absence of an electromagnetic field, Charge Conjugation invariance 
implies that N+ Tq (t) = Ny To (r) and ((iV+) 2 ) = ((iV+) 2 ) for all V. Discrete 
symmetries will be addressed in a future publication. 

We have described, both here and in previous publications [16, 14], the 
utility of a finite-time particle interpretation. In [6] we demonstrated with 
examples how it becomes possible to track the particle production process, 
finding not only how many particles are created, but also when they are 
created, and how they behave after their creation. Having now defined 
finite-volume operators which have well-defined fluctuations, we can also 
say (with controllable precision) where the particles arc created. In a fu- 
ture publication we will show that this makes it possible to treat 'particle 
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creation' and 'vacuum polarisation' effects within the same framework ([8] 
presents a similar treatment for electrostatic fields). In particular, for an 
electromagnetic potential step, modes with evanescent contributions give 
rise to a nonzero charge distribution in the vicinity of the barrier (on a 
length scale A c = ^), providing a vacuum polarisation which partially 
screens the barrier. Meanwhile, the 'Klein modes', present only for a po- 
tential step larger than 2mc 2 , correspond to the standard 'particle creation' 
effect, with created plane wave states (persisting to spatial infinity). Con- 
ventional 'tunneling' methods, based on calculating transmission/reflection 
coefficients, describe only those modes which persist as plane waves either 
side of the barrier. Methods based only on N = J2i{ a t a i + b\bi}, pick up 
the total number of created particles, but this generally must be divided by 
the infinite spatial volume (2tt) 3 S(0) to obtain an average particle density. 
Again, only those modes contribute which persist to spatial infinity. 

Further details of our approach to fermionic QFT are given in [1, 16]. 
We now study a simple example involving a particle horizon; the well- 
known Unruh effect, concerning a uniformly accelerating observer in flat 
Minkowski spacetimc. 

4. RINDLER SPACE: A SIMPLE HORIZON 

The Unruh effect [31, 32] is one of the most studied and most cited 
examples of particle creation, and demonstrates more than any other the 
fact that the concept of particle is observer-dependent. Useful references 
are [7, 33, 8] and the further references therein. An early treatment of 
fcrmions in Rindler space is [34] (presented again in the same authors text 
[8]), while [35] presents a very thorough review of fermions in Rindler space. 
We now consider this problem using the present formalism. We begin 
by showing that the radar time of a uniformly accelerating observer is 
indeed Rindler time, and that the eigenstates of Hi are stationary states 
of Rindler time. This demonstrates the consistency of this approach with 
standard derivations. We complete the derivation in 1+1 dimensions for 
both massive and massless particles. Although the rest of this derivation is 
quite standard, we have also described the spatial distribution of Rindler 
particles, which has received little attention to date. 
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4.1. Radar Time for Accelerating Observers 

Consider first the 1+1 dimensional case. We have 

^(r) S (t(r),*(r)) = (^^,S^) (57) 





1 
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'k = 0' 
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* z 
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FIG. 2. Hypcrsurfaces of simultaneity of a uniformly accelerating observer. 

Consider a point (t ,z ) to the right of the observer's worldline, as shown 
in figure 2. Clearly, the for this point must satisfy: 

t{T + )-t = z -z{ T +) t -t{T-)=z -z{T-) (58) 

From (57), 

e aT+ =a(z + t ) e- aT ~ =a(z -t ) 

from which we deduce that r= ^log( ^°^*° ). For points (t ,z ) to the left 
of the observers worldline (but still in Region U), the roles of r + and r~ 
are reversed, leaving r unchanged. We can therefor drop the subscripts, 
and write: 



r (x) = i-log(^) (59) 
Za z — t 

which is the Rindler time-coordinate, and covers only region U. The hy- 
persurfaces E To are given by t To (z) = ztanh(aro), as shown in figure 2. 
The radar distance p(x) (which is positive by construction) is given by 
p(x) — l log ( a ^ — ^ = I lp g(° " ) I ^ where u = \Jz 2 — t 2 . In (t,u) coordinates 
the metric takes the familiar form: 
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ds 2 = a 2 u 2 dT 2 — du 2 

Two uniformly accelerating observers, each with different a, will have the 
same hypersurfaces of simultaneity, but will 'tick off' these hypersurfaces 
at different rates. The vector field k 11 is given by: 

d d d 
k = z — — t— or =— — in Rindlcr coordinates 

at oz or 

which is the Killing vector field that is used to define positive/negative 
frequency modes in conventional derivations of the Unruh effect. 




FIG. 3. Hypersurfaces of simultaneity of an observer undergoing uniform accelera- 
tion for a finite period of time (in Region U here) who is otherwise inertial. 

To see the significance of the dotted line in region II, consider the 'Finite 
Acceleration Time' case shown in figure 3. In the limit as the 'Acceleration 
Time' approaches infinity, this case approaches that of a uniformly accel- 
erating observer. In this limit the hypersurfaces of simultaneity (which are 
Cauchy) all approach this dotted line in region II (k^ (x) — > there) , and 
a particle horizon forms. This explains the requirement (essential to many 
derivations of the Unruh Effect [34, 8, 7, 36]) that the 'Rindler modes' must 
be those that are zero throughout region II. Further details of the 'Finite 
Acceleration Time' case can be found in [6]. 
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Consider now the 3+1 dimensional case, with x'fdr) = ( &mh ^ aT \ o,0, cosh ^ aT ) ) , 
In this case equations (58) can be replaced by (dropping subscripts) 

{t{T+)-tf = {z-z{T+)f + \ X ^ 

(where \x ^_\ 2 = x 2 + y 2 ) , which applies to both r+ and t~ . Substitution of 
the expressions for t(r) and z{t) and rearrangement gives: 

r ±-(5i L i) + °< 3+ '>) r±+ i ± i=° < 60 » 

where T± = e ar± are the two roots of (60). From this it follows immediately 
that e 2aT =T + T- = just as in the 1+1 dimensional case, fc^(x) is also 
as before, and we can again define u = \J z 2 — t 2 so that the metric can be 
written as 

ds 2 = a 2 u 2 dr 2 — du 2 — dx 2 — dy 2 

reproducing the Rindler coordinates. 

We have now justified the use of Rindler coordinates for studying uni- 
formly accelerating observers, and the derivation of the conventional Unruh 
effect becomes standard [8, 34, 35]. We reproduce the massive and massless 
cases here, restricting ourselves to 1+1 dimensions for convenience. 

By defining 70 = aujo and 73 = 73 (where 70,73 are any matrices satisfy- 
ing 7q = 1 = —73 and {70,73} = 0) we find that T = j a &z an d T 3 = 0. The 
Dirac equation can be written in the form of (17) as: 



■ dtp i , ._ dtp _ n 

t— + -aa 3 = au(-ia 3 — )+awy"ip (61) 
or 2 du 

The inner product in U is written in these coordinates as: 

/•oo 

(i/)\<j>)= / ^(pdu (62) 
Jo 

The M^- term on the RHS of (61) implies, as expected, that H ev is not 
Hermitian. A simple calculation reveals that W ev — H ev — iaa 3 . Hence (61) 
can be rewritten in terms of H\ as: 



dip ' , , ._ dtp. _ n , i _ ,„„. 
i—=H 1 tp = au(-ia 3 — )+au'y u tp--a(j 3 (63) 

OT OU I 
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It is easy to verify that (19) is indeed satisfied. The cigcnstatcs of Hi are 
the stationary states, justifying again the choice made in previous deriva- 
tions. Since the eigenstates of H\ are stationary, then although the Rindler 
observer disagrees that the Minkowski vacuum is empty, he still agrees that 
the particle content does not change with r. 

4- 1.1. Solutions in Region U 
Consider solutions of the form: 

I TPM) = (f(u)<f>++9(u)cl>-)^(u)e- i " r 

where <j>± are basis spinors, satisfying 

<j)\(j) a = 25\ a a 3 (f)±=±<j)± 7o0±=0T 

The factor of 2 is for convenient comparison with the representation given 
in [8]. These are eigenstates of Hi for all r, with eigenvalue oj. 
Substitution into (63) gives: 



u-^ + h+i- )g{u) = mf{u) 



—i 

u V du 2 



(64) 



from which we find that 



(u—u-^-)f(u) = 
y du du' J K ' 

( M T" u T")#(w) = 
du du 



2 2 / W 1 \S 

m u -(---) 
a 2 



/(«) 



(65) 
(66) 



Equations (65) and (66) can be identified with equation (21.66) of [8], or 
with page 4 of [37]. The only normalisable solution is of the form: 

i i/> u (t,v) = {H ( U k (imu)(f) + + HfJ ^{imu)(t>^}e- lWT 



where hI'(z) are Hankel Functions. These states satisfy [8]: 
j j 16aS(uj-w') 

( V>u>l 1>w') = — =^77 

m(l + e a ) 
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We shall return to these solutions after finding a convenient representa- 
tion for the Minkowski modes. 



We have just found a basis for and H^. Together they span all those 
modes that can be seen by the Rindler observer. They span only half of 
TLmi since they do not include any states defined in region II. However, they 
are sufficient to calculate the Rindler number operator N\,r = P^ — P^, 
from which wc can deduce the number of Rindler particles present in the 
Minkowski vacuum, along with their frequency and spatial distributions. 
To define the Minkowski vacuum, we could use the ordinary plane wave 
basis for Minkowski modes. However, it is more convenient for the massive 
case to use an alternative representation of Minkowski modes [8, 35, 34]. As 
a consistency check, we shall rederive the massless limit using Minkowski 
plane wave states. 

The representation of Minkowski modes used in [8] (often called the 
'Rindler Basis') is of the form: 



where 11 tp u (x) = o\i\)u (— x) is a state defined in Region II, and , tpjj 

are states defined in Regions F and P respectively (see [8, 34] for details). 
These various states are defined only in their particular region, so the Dirac 
operator acting on any of these states does not give zero, but gives a distri- 
bution on the light cone through the origin. The coefficients in the above 
expansion are chosen to cancel these distributions. The u^f(x) then form 
an orthogonal basis for H^, an d the v£f(x) form an orthogonal basis for 
Hm (irrespective of the sign of uo). 

In (67) and (68), the states 11 1^^, 1 V^, F V'i ± \' Pl /'i ± ' ) are n °t normalised. 
By rewriting (67) and (68) in terms of states II ip UJl I ip UJ which are each 
normalised to 2tt8(lu — oj'), and choosing N such that u^f(x) and v^f(x) are 
normalised to 2it5(u) — oj'), equations (67) and (68) become: 



4.2. 



The Minkowski modes 



v^ = N(e^ I ^ + e^ F ^+e^ p ^+z II i; ul ) 



(67) 
(68) 



M 



1 




+ F, P, terms) 




v, 



,M 



( e - 7rw/ C orm + i // C° rm + F. p > terms) 



VT+e 
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from which we can immediately extract the relations: 



«L^^H^=£ (69) 
PL^^) = ^t^ (70) 

lL^-^) = 2 -^t^X (71) 
Vl + e— 

<W = ( ip-u,\v u , ) = (72) 
Vl + e^" 

where w is restricted to w > 0, but u/ can take any sign. Alternatively, 

u M 

these coefficients can be extracted using the observation that , "_ 2xu , + 

V l+e a" 

/./.norm 



a 2 (equal to )/)™ m + F, P, terms ) is the only linear combination 

of w^f and vff that is zero in Region II. Although the inner products in 
(69) - (72) have been evaluated on a hypersurface (from (62)) which is 
not Cauchy in Minkowski space (it covers space according to the Rindlcr 
observer), it is still consistent with any inner product in Minkowski space, 
since the Rindler modes are all zero in region II. 
From (69) - (72) it follows that: 

J-oo 27r 1 + c— l + e— 

which represents a thermal spectrum at Temperature T — 2 ^ ks , as ex- 
pected. 

4.2.1. The Massless Case 

In the massless case equations (64) decouple, and we find two indepen- 
dent solutions for each ui, 



\il> u ,i) = M<™)*~*e~ iuT an d |Va,,2> = 0-(au)^ K -5e- <WT (74) 

Each of these is normalised to (iP u ,<t\iPu',(t') = (27t)<5 (T(T '(5((jj — ui'). The mass- 
less limit of I 7 ?/'™™) wou ld yield one specific linear combination of 
and |V>u>,2)- However, we shall find that an d IVV2) eacn l ea d to a 

thermal spectrum, so that we need not restrict ourselves to this one linear 
combination. 
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The plane wave states 7i M also take a particularly simple form in the 
massless case. A basis of is provided by states of the form: 

+e -i P (t-x) and ^^zpit+x) forp>0 

and a basis for is provided by states of the form: 

(jy+e 1 ^-^ and e lp{t+x) for p>0 
This allows us to write: 

where (V>u,,i|4>+e ip( *~ x) > = / (^"^(aa;)^"^ =I(? e *)^-|r(— + |) 

Jo a a a 

and L = (V^.ilVw) = / °° ^ = (2n)5(0). We have evaluated the inner prod- 
uct on the hypersurface t = 0, and we have used [37] (equation (6) page 1) 
in the last line. By using further properties of the T-function we deduce 
that: 

<, = if* » . ' (75) 
as expected. Similarly, 

Jo 27r l + e~ 

and it is straightforward to show that each of N^ 2 is a l so eacn equal this 
value. 

4.3. Spatial Distribution of Rindler Particles 

Although (73) and (75) each show that the spectrum of Rindler particles 
is thermal, they do not tell us the spatial distribution of these particles. 
This can be deduced from (53) and (54). For this purpose, it is clearest to 
work with the radar- like spatial coordinate x= ^l°g( au )- I n these coordi- 
nates our observer is at \ — 0, and |x| represents the radar distance from 
the observer to the point (r,x) (for any r). In these coordinates the inner 
product (62) takes the form = e ax ip^ <f>dx- The massless states 

of equation (74) are of the form e~%~ <fi±e~ lw ( T ^ x h They are plane wave 
states in these coordinates, with the conformal factor cancelling the 
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e ax in the inner product. Define the particle/ antiparticle densities n ± (r,x) 
such that n ± {T,x)dx represents the total number of particles/antiparticles 
within dx of x- From (53) and (54) this is given by: 



« ± (t,x) : 



where 



|= r n ±0J (r, X ) dw 
Jo 1 + 27r 

n w (r,x) = e ax VI norm (T,x)V" 



(76) 



Both are independent of r. 

For the massless case, ViV'w = e~~ ax , so that n + (r,x) = n _ (r,x) and 
both are spatially uniform in \. However, as mentioned in Section 3.5, the 
interpretation of n (t, x) as representing 'particle density' is only accurate 
when averaged over a distance L sufficient to suppress the fluctuations in 
N^(t). From (55) and (56)), a straightforward calculation yields: 



((^V+) 2 ) = ((7V L -) 2 ) 



where 



du doj' sin 2 ((o> — w')L) e~ 



2 r2 roo . 2( 

(77) 



a 2 L 
4ir 2 



sin 2 (a^L) 
d^ ,, 9 _ 9r9 cosh(7r^)£(^) 



/i 2 a 2 L 2 
dX 



1^1 cosh(7r^) + cosh(7rA) 
— log[cosh(7r^)(y / cosh 2 (7r^) + 1-1)] 

TT\J C0Sh 2 (7T^) + 1 




0.2 - 
0.15 

0.1 
0.05 



FIG. 4 A. 



;in (a/j.L) 



FIG 4 B. cosh(-7r/i)L(/i) as a function 
of ii. This function is '"g^+i) at \x = 
as a function of /n, for 0, but quickly approaches its large fj, 



(bottom curve) - and J- (top curve). limit of 



log(2) 
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Figure 4 B shows cosh(7r^)L(/z) as a function of //, while Figure 4 A shows 
s '° 2 ^^ for L = § (bottom curve) i and ^ (top curve). For L»^ the 
integral in (77) is dominated by contributions from small \i. We can then 
approximate cosh(7r/i)L(^) w log ^^ f x ^ w .198. In this limit: 



((^)»)«aL?^+^«.45JV2 (78) 



For i >> i we can approximate cosh(7r/i).L(/z) w °^ - , so that: 

((N+f)~aL^gt = .5Nt (79) 
In either case, we confirm that: 

other wise. 



• ((—£-) 2 ) ~ j^—^- oc X' w hi c h would be infinite for L^O, but is finite 



• The uncertainty in the number of particles in length L, given by y ((iV^ ) 2 ) , 
is proportional to the square root of the number of particles in that length. 
For this uncertainty to be less than the measurement, we must average 
over a length such that N~£ > 1. This is a reasonable requirement - if we 
want a reliably defined particle density, we must average over a volume 
that contains on average more than 1 particle. 

Now we study the massive case. We can show from properties of the 
Hankel functions [37] that n w (r,x) = ti- u (t,x), s ° that the 'principle of 
detailed balance' [35] applies at every point in Rindler space. The relation 
n u> {t, x) ~ n -oj ( t j X) a l so expresses the fact that the distribution of Rindler 
antiparticles exactly matches the distribution of Rindler particles. 

In Figure 5 we have plotted n u (r, x) / a as a function of ax, for m = a 
and u)= |, a and 4a. With c and h included, a length a Ax = 1 corresponds 
to Masses are measured in units of ^§ and frequencies in units 

of |. Normalisation of the masslcss plane wave states to (2tt)S(ui — co') 
represents a norm of '1 per unit length', as is conventional. For massive 
states this interpretation is valid only as a Cauchy principal value, so that 

L 

Ll-^ N ^( T 'X)dx^^ asi^oo.) 

Figure 5 shows that for negative x the Rindler particles are uniformly 
distributed, but for positive x (i- e - to the observer's right) the number of 
particles decreases rapidly. This fact can be understood [35] by writing 
equation (65) in terms of x as: 
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N„ 




-10 -5 5 

FIG. 5. Nu{t,x)/ci as a function of ax, for m = a, and u>- 
4a (most oscillatory curve). 



10 



(lowest curve), a and 



-§ +m2e2ax/=( " + y )2/ (80) 

which takes the form of a 1-dimensional Schrodingcr equation with poten- 
tial V — m 2 e 2ax . For m — the potential disappears and we obtain the 
plane wave states described above. For a — the potential allows V = in 2 
and the mass gap reappears. For nonzero m,a the potential permits states 
of any frequency, but does not allow them to propagate further to the right 
than x = il°g(m) before being exponentially damped. We also see from 
the form of (80) that a change in mass by a factor e al has the effect of trans- 
lating the potential a distance I to the left. This is illustrated in Figure 
6. 

Figure 7 (A) depicts n(r,%)/a as a function of a\, for m = a/10 (right 
curve), a, and 10a, obtained by integrating (76) numerically over < u> < 2a. 
The factor ensures that the number of particles present having en- 

l+e — 

ergy > 2a is negligible (see Fig. 8(B)). We would like to interpret n(r,x)/a 
as the number density at \i an d ^^jja - as 'the number density at x 

a(l+e a ) 

of particles of frequency w'. This again requires caveats regarding fluctua- 
tions, in the same way as in the massless case. Accordingly, we find that the 
particle number density is uniform to the observer's left and negligible to 
the observer's right, and that the position of this transition is determined 
by the ratio — . 

For m— >0 this transition point goes to oo, reproducing the spatial uni- 
formity of the massless limit. However, for non-zero m and realistic accel- 



FERMIONIC QFT IN GRAVITATIONAL BACKGROUNDS 31 



N„ 
a 



J 




\ 1.5 








^0 . 5 


I 



a c 



FIG. 6. N^{r,x)/a as a function of a\, for w = a, and m = a (right curve)and ae 3 
(left curve). 



erations, the particle density at low \ (which is oc a) becomes small, and 
the transition to a negligible density occurs far to the observer's left. 





FIG. 7 A. N(r,x)/a as a function of a\, FIG 7 B. E(x)/a 2 as a function of a\, 
for m = a/10 (right curve), a, and 10a (left for m = a/10 (right curve), a, and 10a (left 
curve) . curve) . 

Consider now the energy density of these Rindler particles as measured 
by the Rindler observer. It is given by E{\) = E + {\) + E~ (\) where 



E ± {X)- 



' w«i u (r,v) da; 



l + e- 



2tt 



Since n w (r,x) is even in uj and is independent of r, it follows that E + (x) — 
E~(x), and is independent of r as expected. Figure 7 (B) shows ^£r- as a 
function of a\, for m — a/10 (right curve), a, and 10a. This is qualitatively 
the same as Figure 7 (A), as expected. 

It is not surprising that the Rindler observer detects the presence of 
particles and of energy; this is fully consistent with detector models. It 
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is changes in the energy levels of the detector that signal the detection of 
Rindlcr particles. However, we should not identify the 'energy as measured 
by an accelerating observer' with the covariant energy momentum tensor 
{Tfiv), which is zero in the Minkowski vacuum, albeit with non-zero fluc- 
tuations. To understand why there is no conflict between E(x) 7^ and 
(Tpi,) = 0, it is instructive to recall the classical connection between T^ v 
and E. For instance ([38], pg 131) if we contract with a 4-velocity 
Up, defined at at some event x, then T^u v represents the "4- momentum 
per unit of three-dimensional volume, '-^r', as measured in the observer's 
Lorentz frame at event x" . This interpretation necessarily involves division 
by a small 3- volume dV , which must be sufficiently small that the observer's 
Lorentz frame remains a good description of his 'rest frame' (and also that 
be effectively constant across V). This in turn requires the volume to 
have dimensions Ax << c 2 a _1 where a is the observers acceleration ([38], 
pg 169). But the uncertainty principle requires that Ax> — (for the un- 
certainty in the momentum to be less than mc). Hence, the connection 

between T^ v u v and the "4-momentum as measured by an observer" ap- 

3 

plies only when a < < - it breaks down precisely when the Unruh effect 
becomes significant. 

We have not discussed either (T^) or the backreaction equation G^ v = 
87r (!),„) which effectively defines it. These issues are discussed extensively 
in the literature (see for example [7, 23, 33] and references therein). We do 
not add to this discussion, except to emphasise that there is no conflict be- 
tween that literature and the observer-dependence discussed here. Though 
the Hermitian operators representing 'observables' may be non-local and 
observer-dependent, these do not affect the evolution of states, which is 
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still local and causal [1, 16]. There is no confict between seeking observer- 
dependent 'observables' and seeking an observer-independent semiclassical 
evolution equation. 

Figure 8 (A) depicts (r, x) /a as a function of — for m — a and a\ = 
—6,-4,-2,0,1 and 2. We have also shown uj < on this plot, to demon- 
strate that n u {T,x) / a = n-< jl {r,x) / a as claimed earlier. This detailed bal- 
ance also implies that the distribution of Rindler antiparticles exactly 
matches the distribution of Rindler particles. For \ = where the Rindler 
particles are uniformly distributed, n u {T,x)/a is completely flat, so that the 
spectrum of particles is exactly thermal there. At \ = we see slight devi- 
ations from a thermal spectrum, becoming more pronounced as \ increases. 
At x = i the deviations from thermal are such that the low frequency modes 
are greatly suppressed. The factor of (l+e - ^ 1 ) -1 in (76) (plotted in FIG. 8 
(B)) suppresses all modes of frequency u) > a, so that overall there are very 
few particles present at x= ^, in accord with the middle curve of Figures 
7 (A) and 7 (B). 

5. FURTHER EXAMPLES OF RADAR TIME 

The radar time of uniformly accelerating observers in 1+1 or 3+1 dimen- 
sions has been presented here and elsewhere [6, 16, 39, 40]. A larger family 
of observers in 1+1 dimensions, including the 'instant turnaround twin' 
(Langevin observer), and the 'gradual turn-around' twin of Figure 3 have 
been described in detail in [6]. In this Section we will further generalise 
these results. In Section 5.1 we present the radar time and radar distance 
of an arbitrary observer in 1+1 dimensional flat space and investigate some 
of their properties. Section 5.2 extends this to arbitrary observers in an 
arbitrary 1+1 dimensional spacetime. Sections 5.3 and 5.4 treat cosmolog- 
ical applications. We shall present the radar time of a comoving observer 
in an FRW universe of arbitrary scale factor a(t) in 1+1 and in 3+1 di- 
mensions, and examine in more detail some examples, including the Milne 
and deSitter universes (for the latter see also [16]). 

5.1. 1+1 Dimensional Flat Space 

The path of an arbitrary observer is completely characterised by a non- 
decreasing function z^\z-), where z± =t±z are null coordinates, and the 

proper time is given by dr 2 =dz + dz- = rf + dz_ on the curve. Define the 
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functions k±(z±) and t ± (z±) by: 



ldz^\z_) , dz h) (z+) 1 



•+ 



r+(z+)= / k+(u)du t~(z-)= / k_{u)du (82) 



o 

The origin is chosen to coincide with the point r = on the curve. On the 
observer's trajectory = k + dz + and t + = r~ = proper time. However, 

these functions are also defined off the trajectory. For points to the right 
of the observer, r 1 *" are as defined in Section 3, and for points to the left 
of the observer the roles of t* 1 are reversed for later convenience. We can 
write the radar time and radar distance as: 



rz+ rz- 

r(:_ : j:r i ( :_ ) ) = 4 ( / k + (u)du + J k-(u)du 



p(z + ,z_) = ±\T + (z + )-T-(z M-I 



/ k+(u)du— I k-(u)di 
Jo Jq 



To make this construction more concrete, consider the case of an inertial 
observer. In this case z+\z-) — k 2 _z_ where k± = = constant, k- 

is the "fc" of Bondi's 'fc-calculus'. =k±z±, so that the radar time is 



i / 1 — v -, / 1 + v t — vx 

which is the time coordinate of the observer's rest frame, as expected. 

For a uniformly accelerating observer we have z^ (z- ) = Y -az > wnere 
the curve has been translated relative to Section 5 so that it passes through 
the origin. This gives k±(z±) = 1±az± , so that: 

2a 1 — az- 2a x+±—t 
p=±- a \\og{{l + az+){l-az-))\ 



which are the translated versions of (59), as expected. 
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Returning to the general case, we can write the metric in coordinates 
(t,/j) as: 

, 9 dT 2 — dp 2 
ds z = - 



k+(z+)k-(z-) 

To writefc±(2;±) in terms of t± involves inverting the expressions (82) for 
T ±( z ±), which is always possible since t±(z±) are both strictly increasing 
functions (for a future-directed timelike observer). It is easy to verify that 
fc + /j_ = 1 for the inertial observer, and that fc + fc_ =e~ 2ap for the uniformly 
accelerating observer. The time-translation vector field, defined in (16), is 
Jpp for all observers. 

Finally, note that V 2 t = = V 2 p, and §§ = §f and §^ = §§ . Hence, if 
we Wick rotate (t,z)^> (s = it,z) and (r, p) — > (£ = it, p) it follows that £, p 
are harmonic, and £,+ip is an analytic (conformal) function of s + iz. 

5.2. Arbitrary 1+1D Spacetime 

All 1+1 dimensional spacetimes can be written in the form: 

ds 2 = n 2 (z +1 z-)dz+dz- = Q 2 {dt 2 - dz 2 ) 

where z± =t±z as before. Much of the previous subsection carries over, 
with the main difference being the parametrisation of proper time, which 
must now satisfy dr 2 = Q 2 dz + dz- on the curve. 
Define: 



fc_(z_) S n(^(z_) > *_) 1 / <fe + )(2! - ) 



dz- 



k + (z + )^n(z + ,z^( z+ )) 



dz™{z+) n(z + ,z { j\z + )f 



dz+ ~ k-(z h \z+)) 
t + (z + )= / k + (u)du r~(z_) = / k-(u)du 



>z Jz^>(z ) 

We have still chosen the origin of the (z+,Z-) coordinate system to lie on the 
observers trajectory, but we no longer require that it coincide with the point 
of zero proper time. The point r = now occurs at (z + ,Z-) = (zq,z^\zq)). 
This is for later convenience. We still have r(z+,z_) = i(r + (z + ) + r _ (z_)) 
and p(z + ,z_) = ||r + (z + ) +r~ (z_)|, and the metric becomes: 

^ CT^ -^ (83) 
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The time-translation vector field is still J^, and the final comments of 
Section 5.1, about the analyticity of the Wick- rotated spacetime, still apply. 

5.3. FRW Universes in 1+1 Dimensions 

Consider a comoving observer in an FRW universe of arbitrary scale 
factor a(t), in 1+1 dimensions: 

ds 2 = dt 2 - a(t) 2 dz 2 = C( V ) 2 (dr) 2 - dz 2 ) 

where C(t](t)) = a(t) and: 

/"* dt ' 

"®=<» + J W) (84) 

This is a special case of Section 5.2. A comoving observer at z = satisfies 
z+\z-) = z- where z± =r)±z, so that 

k±(z±) = C(z±) and t ± (z±)= / C{u)du (85) 

J no 

We have chosen zq = rjo, so that on the observer's trajectory t + =t~ =t 
(which is the observer's proper time). From (83), the metric can be written 

as 

2 + z -)) 2 / , 2 j 2\ (c, a \ 

d3 = C{z + )C{z-) (dT ] (86) 
In this case we can invert t ± (z±) explicitly, to get 

z±=v(r ± ) (87) 

(with r](u) defined as in (84)). This allows us to rewrite (86) in terms of 
r,p as: 

dS= a(T+)a(T-) ( dT - d P) ( 88 ) 

where —r±p. 

To generate insight into these results, we consider 3 examples. 

5.3.1. deSitter Space 

DeSitter space has a(t) — e xt 1 which gives r}{t) = ^e~ xt and C(rj) = 
Note that rj{t) < for all t, so the domain of r(x) will cover only the causal 
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past of the point (r),z) — (0,0), with particle horizons at z= =pr^ = ^-e xt . 
From (85) or (87) we have t ± (z±) = ^log(— \z±), so that r,p are given 
by: 



r = ^log(A 2 z+z_) = ^log(A V - z 2 )) 



= 2l l0g(e 



- 2At -A 2 z 2 ) 



1 e- xt 
"2A l0g( ^ 



+ Az N 



-Xz' 



1 e~ xt 
2A l0g( ^ 



+ A|z|, 



(89) 
(90) 



t 




FIG. 9. DcSitter space in (t, \z\) coordinates. The grey lines represent ingoing and 

outgoing photon trajectories (the outermost one is the particle horizon z(t) = — ) : 
and the black lines are the observer's hypersurfaces of simultaneity for various values of 

TO- 

The metric in (r,p) coordinates, calculated either from (88) or directly 
from (89) and (90), is ds 2 = coshT 2 (\p){dT 2 — dp 2 ). The time-translation 
vector field is given by: 

t 9 d d , . . 

k=— — Az— or =— m (r,p) coordinates 

ot oz or 

Again this is a timelike Killing vector field on the domain of t(x), and 
is spacelike outside this region. This r coordinate, and the corresponding 
Killing vector field, are the same as those used in Gibbons and Hawkings 
original derivation of the thermal deSitter spectrum in 1977 [41], and as 
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in various other derivations since [42, 43] (see also [7] and the references 
therein). The metric takes a more familiar form [44] if we substitute u = 
tanh(Ap) ^ gj vm g = (1 — \ 2 u 2 )dT 2 — 1 _ d ^ 2u 2 ■ Hypersurfaces of constant r 
have been plotted in Figure 9. 

5.3.2. The Milne Universe 

The 1+1 D Milne universe is given by ds 2 = dt 2 — a^t 2 dz 2 . Application 
of (85) now gives t ± =te ±a ° z , so that: 

t(x) = icosh(aoz) p(x) = |<sinh(aoz) | 

In terms of which the metric is ds 2 = dr 2 — dp 2 . It is convenient to drop 
the absolute value in p and to treat p as a spatial (rather than a radial) 
coordinate. Although the cases t > and t < must be treated separately, a 
similar transformation holds in each case. With the coordinate singularity 
at t = moved to the light cone through the origin, the regions t > and 
t < are revealed to be regions F and P of Rindler space, and r,p are 
the inertial coordinates of the underlying flat space. This is illustrated in 
Figure 10. 
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Since t = is just a coordinate singularity, there is no reason for the 
observer's trajectory to terminate there; the 'completed' geodesic will cover 
all t. Clearly, as this observer charts her radar coordinates (sending and 
receiving light from various spacetime points) it will be readily apparent to 
her that regions II and U also exist, and the domain of (r,p) will certainly 
include these regions. The lines t = ±p which bound the original (t,z) 
coordinate system are of no significance in the (r,p) coordinate system. 

Some approaches to particle creation in the Milne universe (see [7] Sec- 
tion 5.3, and references therein) exploit the similarity with the Rindler case 
to derive non-trivial quantum effects in this spacetime. However, an im- 
portant distinction must be made. Region U of Rindler space, expressed in 
Rindler coordinates, corresponds to spacetime as seen by a uniformly ac- 
celerating observer in flat space; it is this fact that gives physical meaning 
to the Rindler coordinate system. In the Milne universe, however, no ob- 
server would have t as their radar time. A comoving observer at the origin 
of the Milne universe corresponds simply to an inertial observer (stationary 
at the origin) in flat space. It is therefore reassuring that the radar time 
of that observer brings us back to Minkowski space, and to the standard 
inertial vacuum. 

5.3.3. Radiation domination 
In this case a(t) — which gives: 



r(x)=t+-z 2 
ds 2 = ±(l + 



p(x) = \z\V\t 

■ )(dT*-dp>) 



Hypersurfaces of constant radar time are plotted in Figure 1 1 . 

As a result of the singularity at t = 0, the comoving observer's geodesic 
starts at finite proper time (chosen for convenience to be zero). The horizon 

now represents the future light cone p = r (or r = ^\J\) of the origin. 
This is not an acceleration horizon, and (unlike the Rindler case) the time 
translation vector field does not tend to zero there. Instead we find that, 
though the direction of the vector field k^{x) becomes null on the horizon, 
its components become infinite in this limit, so that k 2 -^oo as the horizon 
is approached (i.e. as p— >r). 
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t 



Horizon at Izl = 2 4tJ\ 



I z I 



FIG. 11. Hypersurfaces of constant r in an FRW universe with scale factor a(t) = 
\f~Xt, in (i,|z|) coordinates. The domain of radar time covers only the causal future of 
the origin. 

5.4. 3+1 D Cosmologies 

Consider the metric 



(where /(r) — sin(r),r or sinh(r) for spatial sections which are respectively 
hyperbolic, flat, or closed), and consider a comoving observer at the origin. 
Only the radial null geodesies are relevant for determining r for this 
observer, so that we can define z —r/ztr and the 1+1 D results carry 
over almost unaltered. The are still given by (85) or by inverting (87), 
and the metric is given by: 



ds 2 = dt 2 - a(t) 2 (dr 2 + f(r) 2 {d9 2 - sin 2 (6»)<# 2 )) 




) 



where r(r+,r ) = \{rj{T + ) — t)(t )). The examples of Section 5.3 extend 
almost unaltered to the 3+1 dimensional case. 



6. DISCUSSION 



We have presented a formulation of fermionic quantum field theory in 
electromagnetic and gravitational backgrounds that is analogous to the 
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methods used in multiparticle quantum mechanics. The main difference 
is in the particle interpretation, which requires us to consider the entire 
Dirac Sea, as well as any particles which may be present. This approach 
provides a conceptually transparent approach to the theory and a simple 
derivation of the general S-Matrix clement and expectation value of the 
theory. Moreover, it also leads to a consistent particle interpretation for 
all times and any background, without requiring any 'asymptotic niceness 
conditions' on the 'in' and 'out' states. Other advantages include the ease 
with which unitarity of the S-Matrix follows from conservation of the Dirac 
inner product, insights into quantum anomalies, and the fact that Hermi- 
tian extension provides well-defined second quantized operators without 
requiring a complete set of orthonormal modes. We have used the concept 
of 'radar time' to generalise the particle interpretation to an arbitrarily 
moving observer, providing a definition of particle which depends only on 
the observer's motion and on the background, and not on the choice of co- 
ordinates, the choice of gauge, or the detailed construction of the particle 
detector. 

Ever since the pioneering work of Unruh [31] and Davies [32] in 1975, it 
has been known that the concept of particle differs for different observers. 
However, attempts to systematically assign a choice of particle/antiparticle 
modes uniquely to a given observer have not been successful. Such defi- 
nitions have either relied on the existence of suitable symmetries (Killing 
vectors, conformal symmetry etc) or on an arbitrary choice of foliation of 
spacetime into 'space' and 'time' (see for instance [7, 33] and references 
therein). A notable exception is reference [45] which treated arbitrary ob- 
servers in flat space, using a foliation defined uniquely in terms of the 
motion of the observer. However, as well as being only applicable in flat 
space, their foliation is often multivalued, and omits portions of the ob- 
servers causal envelope if discrete changes in velocity are allowed [14]. 

Though particle detector models provide a useful operational particle 
concept, it is inherently circular for a particle detector to be anything that 
detects particles, and a particle to be 'anything detected by a particle de- 
tector'. It is also known [4, 5] that in electromagnetic backgrounds (and 
in cases where particles could already be defined independtly of detector 
models) the predictions of particle detector models are not always propor- 
tional to the number of particles present, even when the detector is inertial. 
We might conclude that a particle concept is only an approximate notion, 
or even that 'particles do not exist' [46]. However, the observer-dependent 
particle interpretation presented here (and in [1]) averts this pessimistic 
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conclusion, and provides a concrete answer to to the question "what do 
particle detectors detect?" 

We hope that the present work, along with [1, 16], has shown the compu- 
tational and conceptual value of working with a concrete representation of 
the Dirac Sea. We strongly support Jackiw's claim [17] that "physical con- 
sequences can be drawn from Dirac's construction" . Care has been taken 
to ensure that the dynamics is kept separate from the kinematics. The evo- 
lution equation is explicitly local and causal. Though the categorisation of 
states in terms of their particle content requires an observer-dependent foli- 
ation of spacetime, this in no way affects the evolution of these states. This 
goes some way towards showing that the foliation dependence of quantum 
mechanics need not conflict with the coordinate covariance of general rela- 
tivity, provided one remembers the important role played in both theories 
by the observer. 

APPENDIX: GREEN FUNCTIONS AND PROJECTION 

OPERATORS 

In equations (39) - (40) the vacuum expectation value of a general 'one- 
particle' operator Aphys (the physical extension of some operator A\ : TL — > 
H) was expressed in terms of traces of projection operators P^ (r) and P^. 
It is common in many textbooks [7, 33, 8] to express expectation values in 
terms of traces of (2-point) Greens functions. This appendix summarises 
the connections between projection operators and Greens functions. We 
also make connection with the "first order density matrix" or "Dirac den- 
sity matrix" of multiparticle quantum mechanics (see [30] pgs 8-10 for in- 
stance), emphasising the role of the negative-energy Wightman function as 
the Dirac density matrix of the Dirac Sea. The observer-dependent par- 
ticle interpretation presented in Section 3, along with the corresponding 
time-dependent vacuum, allows the definition of a time-dependent family 
of 'vacuum Greens functions'. By presenting these as special cases of gen- 
eral state-dependent Greens functions, further clarification can be made of 
the connection between the negative-energy Wightman function and the 
Dirac density matrix. 

In the quantum mechanics of non-relativistic fermions, a common tool 
(see e.g. [30], pg 8-10) is the first order density matrix 7(7-' ,a' ,r, a), defined 
from the normalised many-body wave function $(n,CTi,...,rjv,o'iv) by: 
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l<s>{r',a',r,a) = N ^ J ^ r 2 • ■■dr N $(r,a,.. .,r N ,a N )^*(r',a', ...,r N ,a N ) 

When the wavcfunction is a Slater determinant of one-particle states 4>i, 
this takes the simpler form: 

N 

7$ (r', a', r, a) = ^ & (r , o)4>* (r' , a) 

i=l 

and is referred to as the Dirac density matrix. Consider now the relativistic 
case, and take the multiparticle wavcfunction to be the 'evolved vacuum' 
|vac To (T)) from (22). Then the Dirac density matrix 'at time r' is: 

7|vac T() (t)> {y\r,x\ T ) = { x \r)Vi,r (v\t) 

i 

(which is a 4x 4 matrix, implicitly containing the spin dependence). From 
(37) this is simply the negative energy projection operator P~ (r) in coor- 
dinate representation. More precisely: 

Pr ( T )i , ( x \r) = e(y)7| vaCTo(r)) (y| T ,a;| T )7''( 2 /)V'(y| T )dS A1 (2/) 

The connection between P~(r), 7| V ac T0 (r)> {y, x ) an d the negative energy 
Wightman function S~(x,y) = (0\^(y)tp(x)\0} also follows. To see this, in- 
troduce the Schrodinger picture field operator i/)(x\j:) = ^ Vi^liOViCxIs); 
as in multiparticle quantum mechanics. The connection with Canonical 
methods is straightforward [1]. The field operator can be written in the 
Hcisenberg picture as i>(x) = J2itpi(x)i,/, i ( x \ s ), where S r is an arbitrary 
fixed 'reference hypersurface' on which the Hcisenberg picture states are 
defined. If we choose the 'vacuum' |0) to be |vac To (£ r )) (which is the 
Heisenberg picture representative of |vac To (r))) then we have 

i 

where the subscript r signifies that S~ (x,y) depends on t through the 
choice of vacuum state |vac T0 (r)). Therefore 7| V ac T0 (r)) {y\T,x\ T ) = S~ (x\ T ,y T ). 
That is, the negative energy Wightman function is the Dirac density func- 
tion of the Dirac Sea, and is the kernel of P~(r). The positive energy 
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Wightman function S+ o (x,y) can similarly be written as: 

^n ) ( a; 'y)=) vac ^( s r)|V'(x)^(y)|vac To (i; r ))==^u i , To (x)u i;To (y) 

i 

which is the kernel of P+ (r) . 

In terms of S^ (x,y), other commonly used 2-point functions [33] are the: 

Feynman propogator iSp tT0 (x,y) 
Hadamard function S$ (x, y) 
full propogator iS(x,y) 
retarded propogator Sn(x,y) 
advanced propogator SA(x,y) 

S^l (x\ T ,y\ T ) is the kernel of the 'first quantized' number operator Ni, To (r), 
while iS(x\ T ,y\ T >) is the kernel of the 'first quantized' evolution operator 
U\{t,t'). S(x,y), Sf{(x,y) and Sa{x,v) are independent of t , since they 
do not depend on any choice of state. Furthermore they do not depend on 
the decomposition of solutions into positive/negative frequency modes. 

The Wightman function S~ (x,y) can be generalised to an arbitrary state 
|F(E r )) (which is the Heisenberg picture representative of \F(t))) as: 



= 0(t x - T y )S+ {x,y)~ 9{t v - t x )S To (x,y) 
= S? (x,y)-S-{x,y) 

= Sr ( x ,y)+Sr (x,y) 

= -0(T x -T y )S(x,y) 
= 0(T y -T x )S(x,y) 



S^ F) (x,y) = (F(i: r )\i;(yW(x)\F(X r )) 
= ^2i)i(x)7p(y)=^ lF) (y,x) 



Again, the negative energy Wightman function is just the Dirac density 
matrix of the state in question, and is the kernel of the operator Pyp^ (t) 
defined in Section 3.3. Here S^(x,y) is the kernel of P^(t). 
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